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ABSTRACT 



f 



We present the results of the most sensitive and comprehensive survey yet undertaken 
for radio pulsars and fast transients in the Andromeda galaxy (M31) and its satellites, using 
the Westerbork Synthesis Radio Telescope (WSRT) at a central frequency of 328 MHz. We 
used the WSRT in a special configuration called 8gr8 (eight-grate) mode, which provides a 
large instantaneous field-of-view, about 5 square degrees per pointing, with good sensitivity, 
long dwell times (up to 8 hours per pointing), and good spatial resolution (a few arc minutes) 
for locating sources. We have searched for both periodicities and single pulses in our data, 
aiming to detect bright, persistent radio pulsars and rotating radio transients (RRATs) of either 
Galactic or extragalactic origin. Our searches did not reveal any confirmed periodic signals 
or bright single bursts from (potentially) cosmological distances. However, we do report the 
detection of several single pulse events, some repeating at the same dispersion measure, which 
could potentially originate from neutron stars in M3 1 . One in particular was seen multiple 
times, including a burst of six pulses in 2000 seconds, at a dispersion measure of 54.7 pc 
cm~^, which potentially places the origin of this source outside of our Galaxy. Our results 
are compared to a range of hypothetical populations of pulsars and RRATs in M3 1 and allow 
us to constrain the luminosity function of pulsars in M31. They also show that, unless the 
pulsar population in M3 1 is much dimmer than in our Galaxy, there is no need to invoke any 
violation of the inverse square law of the distance for pulsar fluxes. 
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1 INTRODUCTION 

Extragalactic radio pulsars, i.e. pulsars located in galaxies beyond 
our own Galaxy and its satellites, have yet to be detected. This is 
primarily because they are at least 100 times more distant than the 
majority of the known population of radio pulsars in our Galaxy. 
To date, the most distant pulsars detected have been found in- 
side the Large Magellanic Cloud (LMC) and the Small Magellanic 
Clo ud (SMC), located at d istances of 49 and 57 kpc respectively 
(e.g. lManchester et al.l200 61. Several attempts to detect extragalac- 
tic rad io pulsars have been made in the past. iLinscott & ErkesI 
( Il980() reported the detection of highly dispersed radio pulses from 
M87 with a d uration of 50 ms. The s e findings howeve r were not 



confirmed bv iHankins et al.l ( Il98lh . iMcCulloch eTal] 11981,) or 



contirmea bv IMan 
lTavloretal.l l ll98ll) 
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The brightest rad io bursts known ar e the giant pulses seen 
from the Crab pulsar jCordes et alj|2004l) and a handful of other 
neutron stars. These pulses are very narrow (in some cases with 
structure still unresolve d at 2-ns time resolut ion) and very bright 
(> 10^ Jy) as reported bv lHankins et alj ( l2003l) . Due to their bright- 
ness, giant pulses may in principle be detected in other galax- 
ies and McLaughlin & Cordes (2003) performed searches for giant 
pulses from M33, the LMC, NGC253, NGC300, Fornax, NGC6300 
and NGC7793, b ut with no confirmed detections. More recently, 
Ishat et akl ( I2OIII) undertook a search for pulsars in M33 at 1400 
MHz, also with null results. 

Generally, pulsar surveys use two complementary search tech- 
niques: Fourier-based periodicity searches as well as searches for 
dispersed single pulses. Periodicity searches are poorly optimized 
for sources that exhibit pulses with a wide range of intensiti es, as 
comprehensively discussed by Cordes & McLaughlin' ( l2003h . Sin- 
gle pulse searches can be more powerful than periodicity searches 
if the source emits bright but rare pulses and/or shows a shallow 
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power-law pulse energy distribution. In these searches, individual 
bright, dispersed pulses are recorded and subsequently analyzed to 
establish a possible underlying periodicity, which does not neces- 
sarily show up in a periodicity search. 

In recent years, single pulse search methods have become 
standard practice and are now employed in most of the on-going 
pulsar surveys. This has resulted in two important discoveries. 
The first is the popula tion of rotating radio transients (RRATs, 
iMcLaughlin et alj|2006l) , objects which are characterized by their 
sporadic, bright radio pulses, and w hich possibly represent a subset 
of the total radio pulsar population. IMcLaughlin et ?A\ i2006h show 
that RRATs are among the most luminous radio sources known 
when they emit (L1400 ~ 1-3.6 Jy kpc^), though they are still 
weaker than the giant pulses seen from, e.g., the Crab pulsar and 
PSR BI937+21. RRATs can account for a considerable fraction of 
the active radio-emitting neutron stars in our Galaxy and their high 
luminosity may make it possible to detect them in nearby galax- 
ies. The second discovery is the detection of what appears to be 
a sin gle radio pulse of extragalactic origin: the "Lorimer Burst" 
jLori mer et al. 2007), an extremely bright (5 1400 ~ 30 Jy), ~ 5-ms- 
wide pulse from a high Galactic latitude and with a high dispersion 
measure (hereafter DM) of about 375 pc cm"^, both which sug- 
gest an origin well beyond our galaxy (rf > 1 Gpc). Such a burst 
co uld, e.g., be a ssociated with the coalescence of two neutron stars 



(iHansen & Lvutik ov 2001) or the eva poration of a black hole CReesi 
Il977h . Other detections of this kind JKeane et al]|201lh with other 
observatories and, especially, where the source position can be bet- 
ter localized, would help to elucidate the true nature of this phe- 
nomenon. 

Searching for extragalactic pulsars is motivated by gaining 
a better understanding of galactic pulsar populations, especially 
when the star formation rate and stellar evolution history differ 
from those in our Galaxy. The sensitivity of current radio tele- 
scopes is such that we are in a position to detect only the most 
luminous pulsars in the nearby galaxies of our Local Group. The 
direct comparison of a pulsar location with that of the known su- 
pernova remnants in M3I such as those found by iGelfand et al.l 
(2005), would be difficult due to the small angular size of these 
remnants. However, if the pulsar age is also determined, associ- 
ations may be possible if the supernova remnants are still bright 
jNaravan & Schaudtll988h . This correlation can help determine the 
formation ratio of rotation-powered pulsars versus other manifesta- 
tions of neutron stars like magnetars, and other compact objects like 
black holes in nearby galaxies. By measuring the DM of sources 
in M3 1 and comparing them with the models of free electrons for 
our Galaxy, we can estimate limits on the contributions to the DM 
due to the interstellar medium in M31, assuming that we know the 
contribution due to our own galaxy, with i n the limits of accuracy 
of the NE2001 model of lCordes & Lazid j2002l) (the contribution 
of the intergalactic medium is likely to be negligible). Finally, by 
measuring the rotation measure (RM) it might also be possible to 
obtain an estimate of the magnetic field along the line-of-sight to 
Andromeda. 

The detection of pulsars in other Local Group galaxies is chal- 
lenging, however. In the case of M31, the disk of the galaxy covers 
a significant area of sky ( 190' x 60') and the distance, ~ 778 kpc 
jKarachentsev et alj[2004l) . means that any pulsar signal is likely 
to be extremely faint. An effective search also requires long dwell 
times in order to catch rare, bright bursts. A deep search of M31 
can be achieved by using the WSRT in a special mode called 8gr8, 
which makes optimal use of the array's linea r configuration, in 
combination with the Pulsar Machine (PuMa; IVoute etn]|2002l) 
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Figure 1. A Sloan Digital Sky Survey image of M3 1 . The circles represent 
the half power beam width of the individual WSRT dishes at a frequency 
of 328 MHz, which cover approximately 5 square degrees each. The solid 
lines represent the two main pointings of this survey, PNTOl (south) and 
PNT02 (north), while the dotted line represents PNT03, a complementary 
pointing centered on the core of M31 (see Table[T]for reference). M32 and 
Ml 10 also fit within the beams (see labels on the figure). 



Table 1. WSRT observations of M3 1 . Wp is the number of pointings at each 
position. The start times of the individual pointing are given in Table 7 of 
Appendix A. 



Pointing 


R.A.(J2000.0) 


DEC.(J2000.0) 




Dwell 


Name 


(h.m.s.) 


(d.m.s.) 




Time (hr) 


M31 PNTOl 


00 41 30 


+41 00 00 


4 


8 


M31 PNT02 


00 44 00 


+41 30 00 


4 


8 


M31 PNT03 


00 42 45 


+41 15 00 


2 


8 



backend. The 8gr8 mode provides a ~2.5-deg-wide field-of-view 
(FoV), that of the primary beam of the individual 25-m WSRT 
dishes, as well as the full sensitivity and angular resolution of all 
the dishes combined in a tied-array. This enables the possibility 
of locating new sources with good accuracy (a few arcminutes) if 
the source emits bright periodic pulses or, at least, multiple bright 
bursts at different hour angles. 

In this paper, we present a search of the Andromeda galaxy 
(M31) and its satellites (M32, MHO and a few other dwarf 
galaxies). In §2 we describe the observations and the 8gr8 tech- 
nique; in §3 we present the data analysis and describe the search 
methodology; in §4 we present the search results; and in §5 we 
compare our results to our current understanding of the luminosity 
distributions of pulsars and RRATs in our own Galaxy. Lastly, our 
conclusions are presented in §6. 
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2 WSRT 8GR8 OBSERVATIONS 

We observed two main pointings along the disk of M3 1 (solid cir- 
cles of Fig.[T]l and later also one complementary pointing centered 
on the core of the galaxy (dashed circle of Fig. [TJ. Each of the 
main pointings was observed for 32 h total, divided in to four 8- 
h observations, separated by intervals of approximately one day. 
The complementary pointing was observed for 16 h in two 8-h ses- 
sions at a later epoch. We recorded two polarizations (summed in 
quadrature), each with 10 MHz of bandwid th at a center frequ ency 
of 328 MHz. The digital filter-bank PuMa ("V oute et al.ll2002h was 
used to write spectra with 256 frequency channels and U = 409.6 //s 
sampling, producing 2-* samples per observation. Table 1 summa- 
rizes these observations, which were perf ormed in 2005 Octo ber 
using the WSRT in 8gr8 mode (see also Ijanssen et aLl l2009l) . In 
this mode, 12 dishes of 25 m in diameter, and equally spaced by 
144 m, are used to form an array with the equivalent collecting area 
of a single 74-m dish. 

When the signals of the 12 dishes are coherently added, hav- 
ing corrected for instrumental and geometrical phase delays, we 
obtain an elliptically shaped fan beam. The major axis of the fan 
beam spans the entire primary beam of the 25-m dishes, and the 
minor axis is inversely proportional to the separation of the furthest 
dishes. Due to the regular spacing of the dishes along the East- 
West axis, they produce a grating response on the sky with parallel 
fan beams equally spaced across the primary beam and separated 
by c/(Bv) radians, where B is the projected baseline between the 
dishes, v is the observing frequency and c is the speed of light. 
Each such collection of fan beams are referred to as a "grating 
group" from now on. Since the WSRT has eight independent signal 
chains, which can each have their own separate phase-tracking cen- 
ter, we can simultaneously form and record eight grating groups on 
the sky. 

These grating groups allow us to almost fully tessellate the 
entire primary beam. At a later stage, the eight grating groups can 
be linearly combined in software into A',otai elliptical sub-beams to 
cover the entire FoV of the primary beam by making use of the fact 
that the grating groups rotate relative to the sky during the obser- 
vation. Each of the sub-beams covers an angular size of approxi- 
mately ~ 1.8 arcmin, at an observing frequency of 328 MHz, and 
is spaced in such a way so as to overlap with adjacent sub-beams 
at the half power point. The total number of sub-beams Ntata varies 
per observation between 900 and 1100, and depends on the posi- 
tion of the main beam on the sky (hour angle) and the length of 
the observation. Each sub-beam has a comparable sensitivity to a 
coherent combination of all 12 dishes. Each of these sub-beams is 
searched for dispersed, pulsed signals and single bursts. The sub- 
beams allow us to localize sources with modest precision when we 
have detections from observations made at different hour angles, 
since the different orientations of the fan beams allow one to local- 
ize the source at the crossing points of these beams. Having a large 
number of sub-beams also permits us to distinguish spurious detec- 
tions from real ones because genuine single bursts will be detected 
in multiple sub-beams. 

In the case of periodicity searches, a candidate should appear 
in multiple detection regions in an observation, as shown in the 
upper panel of Fig. |2] The number of beams in which the candi- 
date is detected depends on the intensity of the source. For sin- 
gle pulse searches, localization is more complicated because the 
source will appear in more sub-beams. A genuine astrophysical 
detection should appear in a number of sub-beams A'b such that 

A'max > Nb > A?mi„, whcrC A'raax ~ l/8Motal and A'„„n ~ l%A'total 
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Figure 2. An example of how sources can be localized using the 8gr8 
method. The panels show two examples of beam plots. The upper 
panel shows the original and confirmation observations in which PSR 
J1937+2950 was detected using a periodic ity s earch (position shown with 
a cross, figure taken from Janssen et al.l 120091) . These observations were 
each 2 h long and were done at hour angles of HA| = 138.81° and 
HA2 = 201.31°. In the lower plot we show the original and confirmation 
observations in which single pulses from the known pulsai' PSR B 1946+35 
were detected (position shown with large cross). These observations were 
done at hour angles HAi = 174.45° and HAt = 244.45° respectively. In 
both panels solid dots represent sub-beams in which the pulsar was detected 
in the first observation while the open circles represent detections made in 
the confirmation observation. The sizes of all the circles are proportional to 
the S/N of the detections. In both cases, the two observations at different 
hour angles allow us to accurately locate the source. 



with the sub-beams making a specific pattern within the primary 
beam as shown in the lower panel of Fig. [2] Here we show the de- 
tection of 8 pulses in total, divided into two groups of 4 pulses, each 
group from a different observation performed at a significantly dif- 
ferent hour angle. A detection appearing in too many sub-beams or 
too many fan beams (A^b > A^max) is likely to be impulsive radio 
frequency interference (RFI). In contrast, a detection that appears 
in too few sub-beams (A'b < A'min) is likely to be a statistical fluctu- 
ation. 
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3 DATA REDUCTION 

As we describe in detail below, we searched for both periodic and 
sporadic signals. We first discuss the common elements of the anal- 
ysis. 

The NE2001 model jCordes & Laziol2002l) predict.^] a DM = 
55 ± 10 pc cm"^^ for a pulsar at the outer edge of our Galaxy, along 
the line-of-sight to M31, which lies well out of the Galactic Plane 
{b = -21°). To detect sources located well inside M31, we searched 
a broad range of trial DMs from to 350 pc cm"'. This DM range 
covers sources located within our Galaxy (DM < 45 pc cm"'), as 
well as objects that might be located in M31. Our upper limit for 
the DM trials was set due to the fact that at this high DM the radio 
pulses are likely severely scattered at 328 MHz. We consider that 
any value of DM larger than ~ 45 pc cm"' may indicate a source 
location in M31; assuming that the electron distribution in M31 is 
not too dissimilar to the one in our Galaxy, this DM range would 
allow us to detect sources reasonably deep into M31, even though 
it is i nclined at only about 13° to our line-of-sight jSimien et al.l 
Il978l) . For example, an excess in DM of 10 pc cm"' (above that 
contributed by our Galaxy) would imply a location about 1 kpc 
inside the outer edge of Andromeda. 

Due to the large data volume, each 8-h observation was sub- 
divided into chunks of 1 h for further processing. For each 1-h 
chunk, we prepared partially dedispersed time series for each of the 
8 grating groups. This was done using a modified version of the tree 
dedispersion algorithm developed by Tavlor ( 1974), which divides 
the bandwidth of each fan beam into smaller frequency sub-bands 
where we assume the dispersive delay to be linear. From these, we 
form time series for all trial DMs (see below), each consisting of 2^' 
samples. Lastly, a geometrical correction was applied before com- 
bining the data sets from each of the 8 fan beams into a collection 
of sub-beams spanning the whole hour. 

Each of these sub-beams was searched for periodicities and 
bright pulses. We searched 549 trial DMs, with step sizes of 0.348, 
0.697 and 1.392 pc cm"' respectively, each corresponding to the 
spacing below the first and second "diagonal DM'0, located at 
DM = 83 and 249 pc cm"' respectively. The DM step-sizes were 
chosen so that the maximum error in the delay between the high- 
est and lowest frequency channels was equal to the sampling time. 
This search used the PSE code from the Parkes High-Latitude Sur- 
vey tedwards et ai]r2001.) , with single-pulse search extensions de- 
veloped by us (see §3.2). 

3.1 Periodicity search analysis 

We performed a periodicity search following the same procedure 
used by Jans se n et al. I i2009l) , but adapted to the specific proper- 
ties of our data. We calculated the power spectra for each DM trial 
of each sub-beam using a Fast Fourier Transform (EFT), remov- 
ing frequencies known to be related to REL Interpolation of the 
spectra was used to recover spectral features lying between Fourier 
bins (Ransom et al. 2002), repeating the process after summing 2, 
4, 8 and 16 harmonics. All the spectral features with a S/N > 7 
were recorded, later keeping only the signals that were detected at 
multiple DMs and appearing to be harmonically related. This was 

' The systematic enor on this predicted maximum DM is ill-defined, but 
from other examples it is conceivable that this value could be off by 20% or 
more. 

^ At the diagonal DM, the dispersion delay across a frequency channel is 
equal to twice the sampling time. 



done to find the most significant candidates, and to determine their 
DM and fundamental frequency. For each candidate, the initial pe- 
riod from the spectral search was optimized by folding the data at a 
range of periods, assuming a maximum error in the pulse frequency 
equal to the width of one Fourier bin. No acceleration searches were 
done. Typically there were a few tens of candidates per sub-beam, 
i.e. about 30,000 candidates per pointing. Those combinations of 
DM and spin frequency that produced the highest S/N candidates 
were stored for later inspection. This process was performed on 
each of the Motai sub-beams. Once this analysis was completed, a 
list of candidates from the sub-beams was collated and those can- 
didates with high S/N appearing in more than one sub-beam were 
stored. Later these detections were compared with the lists from the 
other follow-up pointings, and those candidates that matched with 
similar periods, DMs and S/N were inspected by eye. 

3.2 Single pulse search analysis 

We searched for single pulses following two selection criteria: i) 
single, relatively bright pulses (S/N > 7.0) or ii) multiple, fainter 
pulses below this threshold, but repeating at the same DM. We first 
discuss the general algorithm to search for single pulses, and then 
discuss the details of the two selection criteria. 

As mentioned in §3, for each sub-beam dedispersed time se- 
ries are generated at all trial DMs. Finding single pulse candidates 
simply consists of recording those individual events that show an 
amplitude above some threshold of S/N, which is statistically sig- 
nificant and reduces the number of candidates due to RFI. An ex- 
te nsive description of the meth od for these searches can be found 
in lCordes & McLaughliij|2003l . 

We developed a code to extract individual pulses from each of 
the time series. We first determine the statistics of the noise in the 
time series, which each contain Wsamp samples. This is done by cal- 
culating the mean and the rms using an iterative process where the 
strongest pulses are removed until one reaches the expected num- 
ber of spurious detections (A^noisc) for a Gaussian distribution with 
probability P{x), given by A'noisc = N^.^n-ipPi> 5<r), and then search- 
ing again for pulses above a 5cr threshold. This process keeps the 
mean and the rms from being biased towards very strong individ- 
ual pulses. The S/N of a detection in the time series depends on the 
intrinsic width of the pulse W as S/N oc H'"''^, where the pulses 
with the highest peak S/N will be those with a width similar to the 
time resolution of the timeseries. Therefore we explored a range of 
widths, such that Wj > 4, by adding between 2" and 2^ adjacent 
samples, corresponding to a range in time resolution of 0.4096 ~ 
52.4288 ms. Each of these smoothed time series was then searched 
for single bright pulses, recording all the pulses that showed a peak 
S/N > 5. For each pulse detected we record the DM, the arrival 
time of the pulse relative to the start of the observation and its S/N. 
Hereafter we will call this information DM-S/N-time arrays. The 
S/N threshold was chosen based on the expected number of false 
detections A'(< threshold) above that minimum S/N by assuming 
that the noise is purely Gaussian. For a time series the number of 
expected false detectio ns, above a given S/N, can be c alculated us- 
ing a relation shown bv lCordes & McLaughlinI ( l2003h : 

A'(> threshold) a 2A'samp exp (-^^/2) dx. (1) 

For a full-resolution, 1-h data set, we expect A'(> 5cr) ~2700 false 
detections from all the time series required to cover the 549 DM 
trials. 
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Figure 3. Diagnostic plot from the single-pulse search, for a single sub-beam and a DM range - 83 pc cm" . The top panel shows a DM-S/N-time array, 
with a time resolution of 3.2768 ms, which corresponds to the addition of 2^ raw samples of 409.6 fis. The data ends before t=3600 s because we analyzed 
a timeseries which is the closest power of two number of samples, which for this case is 3436 s long. The sizes of the ckcles are proportional to the S/N, 
with the highest S/N ~ 7.9 associated with RFI at low DM. Those detections that matched our RFI filter have been excised from the figure. The second panel 
from the top shows the time series collapsed along DM in order to more clearly show at which times strong pulses occun'ed. The lower left panels show two 
histograms. The top one is the time series collapsed in time, clearly showing a peak at DM ~ 54 pc cm"^. The lower one shows a histogram where time is 
collapsed to show the number of pulses versus DM. The horizontal lines in the bottom left panels show the mean and the 3 and 6-it limits of each distribution. 
The histogram in the bottom middle pail of the figure shows the distribution of pulses with S/N > 5. The dashed line in this histogram shows the expected 
noise distribution. This same pointing position was observed for another 3 1 h without recording the same burst profusion. 



After searching the dedispersed time series for pulses above 
the S/N threshold, we created diagnostic plots like the one shown 
in Fig.[3]to better assess whether any of the signals were astronom- 
ical in origin. Such a plot was made for each of the sub-beams. 
When sorting/inspecting the diagnostic plots, we searched for both 
i) single bright pulses (S/N > 7.0) and ii) multiple, weaker pulses 
recurring at the same DM. As can be seen in Fig|3] the signature 
of a single bright pulse appears in the top panel of the figure as a 
large circle (e.g., one can see RFI between t = 500 and t = 1000 
s), while the signature of multiple weak pulses at the same DM ap- 
pears as a peak (in this case at DM = 54 pc cm""*) in either of the 
two histograms in the lower left area of the plot. 



3.3 Single bright pulse search 



We recorded all the single pulses with a S/N > 5 in the DM-S/N- 
time arrays for each of the smoothed time series, keeping only the 
information corresponding to those pulses that were detected in a 
reasonable number of beams (l/SNtavA > > l%A'toiai). For each 
candidate, we recorded the DM, S/N, time of the event, the number 
of sub-beams in which the pulse was detected, and in which sub- 
beam the detection was strongest. We then searched for possible 
associations, i.e. other pulses occurring at the same DM at later 
times or in other observations. Those events that had high S/N (> 
7.0) but that did not show any associations were classified as single 
bright pulses. 
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Figure 4. Minimum detectable flux density for single bright pulses as a 
function of the number of samples added, assuming a constant peak ampli- 
tude, so that sensitivity appears to improve with increasing width. The raw 
sampling time of our data is ?s = 409.6 /iS and a minimum thi'eshold of 
S/N=5 was used. Each horizontal line represents the minimum detectable 
flux for each smoothed time series adding from top to bottom, 2" to 2' 
samples. 
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Figure 5. The main panel shows only detections recorded with a DM > 
10 pc cm"' (to avoid low-DM RFI), following the procedures described 
in §3.3. All the pulses presented here appeal' dispersed, were detected in 
many sub-beams, and showed adjacent friends. The inset shows the number 
of expected false detections for purely Gaussian noise, i.e. not including 
detections due to RFI. The tail of this distribution is shown with the dashed 
line in the lower left corner of the main figure. Note the vertical logarithmic 
scale. 



3.4 Multiple pulses at the same DM 

To search for multiple, potentially weaker events occurring at the 
same DM we performed a statistical analysis, searching for any ex- 
cess of S/N in the S/N-DM space. Any excess of S/N due to real 
pulses would create peaks like those shown in the lower left panels 
of Figure [3] Detections showing features of this kind, and appear- 
ing in a reasonable number of sub-beams (A'b > l%A'totai) were 
recorded in a list for further inspection. Since there are 900-1 100 
sub-beams per hour of observation, this produced a very large num- 
ber of candidates. In order to reduce the number of candidates to 
inspect, we filtered the list by searching each candidate for "adja- 
cent friends", i.e. other dete ctions located at adjacent DM trials and 
occurring at the same time. lCordes & McLaughlinl j2003l) quantify 
how the S/N decreases when a pulse is detected at an incorrect DM. 
Narrow pulses will be detected in few DM bins, therefore will show 
few adjacent friends, while broader pulses are detected along many 
DM trials. An example event with adjacent friends is shown in the 
lower left panels of Fig. [S] This criteria was implemented to esti- 
mate the number of adjacent friends a signal of a given DM should 
have, given the S/N of two adjacent detections. Once the data was 
filtered in this way, the remaining diagnostic plots were inspected 
by eye. 



3.6 Sensitivity of WSRT in 8gr8 mode 

For periodicity searches, the minimum detectable flux den- 
sity is given by the mo dified radiometer equation as follows 
jLorimer & Krameill2005l) : 



s/N,„„r,y, / IV 



G^jri^i^ V P 



W 



(2) 



where p is an efficiency factor ~ 1.3 related to the 2-bit sampling. 
S/N,nin is the detection threshold, r^y, is the system temperature in 
K (rsky + Treceiver)! G IS the gain of the telescope (in K Jy" ' ), «p is the 
number of polarizations. A/ is the bandwidth, P is the pulsar period 
and W is the intrinsic width of the pulse. Using Tsys/G = 140 Jy, 
A/ = 10 MHz, /ip = 2, a dwell time of =3600 s, and assuming 
a mean duty cycle of 5%, we can reach a minimum flux density of 
~0.25 mjy for a S/N=7 detection. 

Similarly, for the single pulse searches, our sensitivity is given 
by the radiometer equation, assuming matc hed filtering to a pulse 
that is temporally resolved, as follows dCordes & McLaughlinl 
I2OO3I) : 



GV«pA/iy' 



(3) 



3.5 Radio frequency interference excision 

To excise RFI from the single-pulse diagnostic plots, we inspected 
the DM-S/N-time arrays of each of the 8 grating groups. Events 
aligned in time and appearing along all trial DMs in all 8 data sets 
were flagged as RFI and recorded in a mask file, in order to excise 
those time bins when the single pulse analysis was performed. 



The minimum detectable flux density depends on the pulse width 
as W"''^, assuming a constant peak amplitude, so the sensitivity 
appears to improve with increasing width. As described in §3.2, 
we smoothed the data before recording the peak fluxes in order 
to detect pulses of different widths. The sensitivity curves for each 
smoothing are shown in Fig.|4]for a threshold S/N,^i„=5. It is impor- 
tant to note that the sensitivity decreases for broader pulses because 
we have effectively added more bins above the threshold. 
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Table 2. Summai'y of the pulses shown in Figure |3] These pulses all peak 
in S/N at DM = 54.7 pc cm"' and were detected in the same smoothed 
time series (pulse width of ~ 3.2 ms). The arrival times are relative to MJD 
53707.72175926 (topocentric). 



DM 


S/N 


Arrival Times 


# Sub-Beams 


(pc cm"') 




(s) 




54.7 


5.1 


623.0975 


49 


54.7 


5.5 


658.2706 


13 


54.7 


5.8 


849.0459 


24 


54.7 


5.7 


1345.6613 


121 


54.7 


5.5 


1616.5904 


22 


54.7 


5.8 


2617.5119 


53 



4 RESULTS 

4.1 Periodicity searches 

As described in §3. 1, we performed a periodicity search on each 1- 
h chunk of data. No periodic sources were found using this method. 
We note that there are no known foreground pulsars in the direction 
of our observations. The closest known pulsar is PSR B0053+47, 
located about 7.11° away from the center of our nearest pointing. 
In this search, the minimum detectable period for a DM ~ 50 pc 
cm"' pulsar is roughly a few milliseconds. 



4.2 Single pulse searches 

The single-pulse search produced a number of interesting candi- 
dates. As explained in §3.3 and §3.4 we sifted our single-pulse 
candidates to look for both bright individual pulses (S/N > 7.0) 
and for multiple, weaker pulses occurring at the same DM. We 
discuss the results of each sifting method below. In Figure |5] we 
show a histogram with the distribution of the S/N for the brightest 
pulses reported here. All these pulses were detected according to 
the methodology described in §2. The excess of detections above 
7o" that match our selection criteria, and that are shown in this plot 
implies that these pulses may have an astrophysical origin in the 
absence of RFI or instrument related noise. 



4.2.1 Single bright pulses 

We have detected a number of signals, not obviously associated 
with RFI, with S/N > 7, DM between 10 < DM < 315 pc cm"\ 
and pulse widths of 1 < W < 27 ms, as calculated from the (down- 
sampled) time series they were detected in. We present a table of 
these detections, as well as corresponding diagnostic plots for se- 
lected candidates with high DM, in Appendix A. 



4.2.2 Multiple pulses at the same DM 

Our search for multiple pulses at the same DM revealed a few 
events that satisfy the conditions described in §3.3. The most con- 
vincing detection is shown in Fig. [3] in which 6 weak bursts (5. 1 < 
S/N < 5.8) were recorded during a single hour. The details of each 
of these detections are listed in Tabled This is the most tantalizing 
detection made in our survey; we did not observe any other candi- 
date with the same burst profusion. At this DM, other pulses were 
detected in the remaining 3 1 h of data, but these occurred only spo- 
radically. 

For a Gaussian noise distribution, and using the dwell times 



14 - 
12 - 
10 - 



Theoretical S/N 
Measured S/N 



2 ft 



DM ( pc cm" ) 

Figure 6. Summed S/N of the 6 pulses shown in Table |2] plotted over the 
theoretical S/N curve for a dispersed radio pulse with these properties and 
observing parameters. The eiTor bars in the ,r-axis correspond to half the 
DM spacing while the errors in the j'-axis were calculated from the differ- 
ential form of Eq. (3). 



of our observations, the sampling time and the number of beams, 
we estimate that the probability of having a ~ 5. So" detection, with 
adjacent friends is about < 10"' and that the probability of having 
6 bursts with adjacent friends is much lower (<K 10"'). 

Using the PRESTO tool rrat_period, which performs a 
brute-force period search given a list of pulse arrival times, we 
determined two candidate periods based on the burst times listed 
in Table |2] This analysis assumes that the pulse period is > 50 
ms. If the actual period is less than this, then the time gap be- 
tween the pulses is too large to unambiguously determine a pe- 
riod. Using all 6 arrival times, the best candidate period is 0.23294 
s, implying the following number of rotations between the 6 
arrival times: 150.9977, 818.9964, 2131.9647, 1163.09577, and 
4296.9455. Dropping the first, weakest burst, the best candidate 
period is 0.29578 s, implying the following number of rotations 
between the 5 arrival times: 644.9923, 1679.0073, 915.9843, and 
3384.0160. To check how robust these candidate periods are, we 
performed numerous trials in which we adjusted the pulse arrival 
times to arrive randomly within a 10-second window after the true 
arrival time. We then re-ran rrat.period. Most of these trials re- 
sulted in candidate periods that fit the randomized arrival times just 
as well as for the real data. Thus, we caution that the candidate pe- 
riods we present here do not conclusively indicate that these bursts 
have an underlying periodicity. 

To better test whether these bursts are indeed astronomical in 
nature, we compared the measured S/N of the sum of all 6 pulses, 
as a function of DM, with the expected theoretical curve (following 
Eqs. 12 & 13 given in Cordes & McLaughlin 2003). Figure|6]shows 
that the summed signals match the theoretical curve well. To make 
this summed profile, we added the time series of each pulse, using 
the peak of each burst as the fiducial alignment point. The resulting 
best pulse profile is shown in the top panel of Fig. |7] From this, 
we can conclude that the pulse width is comparable or shorter than 
the sampling time of these time series, i.e. < 3.2 ms. The bottom 
panel of Fig. [7] shows the time series of adjacent DM trials, after 
being corrected for dispersion and spaced in DM as in our original 
search (see §3). The figure clearly shows the dispersed nature of 
these bursts. 

In total we have detected 31 other candidates where weak 



© 2012 RAS, MNRAS OQO.fTlflS] 



8 E. Rubio-Herrera, B. W. Stappers, J.W.T. Hessels & R. Braun 



10 ms 




0.2 0.4 0.6 0.8 

Best Pulse Profile (s) 
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Figure 7. Summed pulse profile of our best candidate. The upper panel 
shows the pulse profile obtained after adding the 6 pulses listed in Table 
|2] The lower panel shows how this summed pulse widens and decreases in 
peak S/N at adjacent DM trials both above and below DM=54.67 pc cm"^ 
(thick line). Due to the very low S/N of the individual pulses, a frequency- 
time waterfall plot showing the dispersion sweep is not informative. 



42.5 




39.5 ' ' ' ' ' ' ' — ' 

12 11.5 11 10.5 10 9.5 9 

RA(deg,J2000) 

Figure 8. Beam plot of the bursts detected at DM = 54.7 pc cm"^^ . The size 
of the circles is proportional to the S/N of the combined detections. The 
highest S/N regions are located in the upper left side of the plot (see crosses) 
and are RA = 11.62°, DEC = +41.44° and RA = 11.19°, DEC = +41.69° 
respectively. The large circle represents the approximate size of the primary 
beam. 

4.2.3 Localization of the candidates 

Due to the sporadic nature of the single pulses, and the linear na- 
ture of the WSRT array, we were not able to unambiguously locate 
any of the sources reported here within the primary beams. To un- 
ambiguously locate a source, we require at least a handful of bright 
single pulse detections at different hour angles, as explained in §2. 
We show an example of this procedure for the pulses of the candi- 
date M31-BC1 (see Table[3ll. Combining only the S/N of the detec- 
tions listed in Table |2] we obtain two best locations for this source 
which are shown in Figure[8] Applying this procedure and includ- 
ding all the detections, we find the best locations for the candidates 
described in §4.2.2 and listed in Table[3] 



pulses repeated with a lower burst rate than the candidate described 
above. From these 31 we report here those with the highest burst 
rate. Their characteristics are shown in Table[3] These sources have 
been detected throughout the area covered by our observations and 
some were found where our observations overlap (see Fig.[TJ. They 
were identified as promising because of a significant overabun- 
dance of pulses at a specific DM. The S/N of these pulses range 
from 5 < S/N < 7. The typical burst rate detected with the WSRT 
for these candidates was between 1-3 per hour. The DMs of these 
detections are generally consistent with that of an object located 
at least in the outskirts of our Galaxy and potentially in M3 1 . For 
each of the sources listed in Table [3] we calculated the most likely 
location within the primary beam by combining the data from the 
individual pulses with the highest S/N in each observation. In Ta- 
ble |BT] of Appendix B we present the brightest multiple pulse de- 
tections we recorded. These events repeated more than once at the 
same DM and with a similar pulse width. 



4.2.4 Follow-up observations 

The candidates shown in Table [3] were followed up first by obser- 
vations using the standard WSRT tied-array mode. These were per- 
formed at the same observing frequency and using the same num- 
ber of frequency channels as in our 8gr8 observations. At this fre- 
quency, the WSRT tied-array beam has a radius of ~ 3', about half 
the size of an 8gr8 sub-beam. Each position was observed for ap- 
proximately 2 h. These observations did not confirm any of the can- 
didates. The low S/N of most of these detections makes differenti- 
ation between real pulses and spurious (noise-related) pulses very 
difficult, even though the pulses appear dispersed and are recorded 
in multiple sub-beams. This might introduce systematic errors in 
the inferred positions that are larger than the 8gr8 sub-beams, and 
hence it is entirely possible that we were simply not pointing in the 
right direction during the confirmation attempts. 

In order to circumvent this problem, we decided to attempt 
confirmations with an instrument that provides larger beam size 
and greater sensitivit y: the Green Bank Telescope (GBT) and Pul- 
sar Spigot backend jKaplan et aLl l2005h at a center frequency of 
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Table 3. The positions of our best four candidates detected using tlie criteria desciibed in §4.2.2 for multiple pulses detected at the same DM. 
For each candidate we hst the right ascension (RA) and declination (DEC) of the best positions (found using the method described in §2), 
the number of bursts detected per hour N (incorporating the full data set), the DM, the inferred distance (d) using the NE2001 Galactic free 
electron model and the width of the individual pulses W corresponding to the time resolution at which the pulses had maximum significance 
in the search. The errors in position for RA can be up to 0.5 arcmin while in DEC they are 1.5 arcmin. The numbers in parentheses represent 
the uncertainty on the least significant digit in the DM. 



Name 


RA(J2000) 


DEC (J2000) 


N 


DM 


d 


W 




h m 


d m 


h-' 


pc cm"^^ 


kpc 


ms 


M31-BC1 


00 43 


+40 22 


2.7 


54(1) 


>2.9 


~ 3.3 




00 44 


+40 39 












00 43 


+40 58 












00 46 


+41 26 












00 44 


+41 44 










M31-BC2 


00 39 


+41 38 


2.7 


56(1) 


> 3.2 


-3.3 




00 39 


+41 14 












00 39 


+40 30 










M31-BC3 


00 37 


+40 24 


3.0 


63(1) 


>48 


~ 3.3 




00 47 


+40 53 












00 40 


+40 40 










M31-BC4 


00 45 


+40 54 


1.0 


147(2) 


>48 


-6.6 




00 43 


+40 45 











43 




39.5 ' ' ' ' ' ' ' ^ 

12.5 12 11.5 11 10.5 10 9.5 9 

RA(deg,J2000.0) 

Figure 9. The medium-sized circles show the area covered by follow-up 
GBT observations. The FoVs of the original WSRT observations are repre- 
sented by the large circles. The black dots show the position of the candi- 
dates listed in Table[3] 

350 MHz. Each of the best positions was observed for 2 h, for a 
total campaign time of 16 h. The pointings observed with the GBT 
are shown in Figure |9] We used a bandwidth A/ = 50 MHz with 
1024 lags and a sampling time of 81.92 yus. The GBT and Pulsar 
Spigot provided roughly a factor of two higher gain and five times 
the recorded bandwidth, meaning an improvement of a factor of 4 
in sensitivity over the original 8gr8 observations. At this frequency 
the beam size of the GBT is .6°. We reduced the se observations 
using the PRESTO package ( Ransom et al. I l2002|) . applying both 
periodicity and single pulse searches. From the burst rates cited in 



Table [3] we were expecting to detect at least one of these candi- 
dates, but no confirmations were made. Assuming that these can- 
didates are indeed from true astronomical sources, the lack of de- 
tections can be attributed to two reasons: i) incorrect localization 
due to the low S/N of the original detections and ii) the possibility 
these sources manifest in rare, sporadic clumps of bursts, requiring 
a longer dwell time to detect. 



5 DISCUSSION 

Our search of M31 can be comp ared with, e.g., the Arecibo 
350-MHz search of M33 by McLau ghlin & CordesI (2003): both 
searches detected some sporadic single pulses with DMs consistent 
with an origin in the target galaxy. To compare the two searches, 
we consider that: i) M33 is ~120 kpc more distant than M31 
dKarachentsev et al] 1 20041) . making any source 1.3x fainter than 
it would be in M31, (ii) the Arecibo radio telescope is 10 times 
more sensitive than the WSRT at these frequencies, and (iii) the low 
mass, different metallicity and star formation rate of M33 would 
likely produce a smaller pulsar population than that of our own 
Galaxy and M3 1 . Taking into account these difference s, it turns out 
that ou r results are similar to those found by McLaughlin & CordesI 
( |2003|) for M33, as we demonstrate below. 

5.1 Detectability of pulsars and RRATs in M31 

Here we consider the detectability of pulsars in M3 1 using the pe- 
riodicity and single pulse search techniques. This, in turn, is used 
to place constraints on the total pulsar population of M3 1 . To do 
this, we constructed a simple population synthesis code to model 
the basic properties of a galactic pulsar and RRAT population. Us- 
ing the sensitivity thresholds listed in §3.6 and our synthesis code, 
we model different populations to find an average expected number 
of detectable pulsars and RRATs. For periodicity searches we used 
integration times of 1 h, while for single pulse searches we used 
observation times of 8 h. 
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5.1.1 Modeling pulsars 

In our pulsar population model, we simulate the distribution of pul- 
sar periods N(P), luminosities as a function of pulse period L(P), 
and the distribution of pulse energies for each pulsar. Each synthe- 
sized population contains about 40,000 pulsars, w hich is based o n 
simulations of our own Galaxy (see, for example, lLorimedl2009h . 
In other words, we assume that the number of pulsars in M31 and 
our Galaxy is similar. We do not simulate millisecond pulsars be- 
cause they do not represent the bulk of the population of pulsars in 
our Galaxy and, in any case, we had poor sensitivit y to them. For 
the pe riod distribution we use the function found by iLorimer et al.l 
l l2006i) : 



N(\og P) = A exp 



(logP-Bf 



2C2 



(4) 



For the distribution of luminosities we derived an empirical relation 
from a sample of 663 pu lsars with known 5400, w hich we obtained 
from the ATNlH catalog I IManchester et al.l2005b . We then fit a log- 
log distribution with a normal distribution given by: 



\og(L) = D\og(P) + E + £:. 



(5) 



Where ^ is a normal distribution centered at the origin with stan- 
dard deviation 0.6. The parameters of our population models are 
summarized in Table |4] An example of the period and luminosity 
distributions for one of our models is shown in Fig. 1101 

For each pulsar, we generated a random pulse-energy distri- 
bution, which was either i) normal, ii) log-normal, or iii) a power 
law. These are th e distributions observed in known pulsars (e.g. 
ICaims et al. I I2OOII and Weltevrede et al. 2006). For power law dis- 
tributions, the index a was assigned randomly, with values between 
I < a < 4. For the normal and log-normal distributions, the mean 
was assumed to be the mean flux density of each pulsar calculated 
from the luminosities given by Eq. (5), using the definition of lumi- 
nosity L328 = Si2>,d-, with d = 778 kpc (Karachentsev et al. 2004, 
the distance to MSI). For the normal distributions, we randomly 
assigned a cr to simulate the distributions of pulses in the range 
100 > (T > 1. For the log-normal distribution, we followed the 
same approach, assigning randomly a a that would produce a range 
of distributions. For each pulsar we also calculated the number of 
pulses that we should observe in one of our 8-h observations (given 
by A'p = t-JP). To build these models, the luminosities of pulsars at 
328 MHz (L328) have been scaled from those at 400 MHz, follow- 
ing the relation S oc v"", where S is the flux at 328 MHz, v is the 
observing frequency an d a is the spectral in dex, assumed to be on 
average a= -1.82 from lMaron et alj ( [2000h . 



5.1.2 Modeling RRATs 

RRATs have been relatively recently identified as a class of radio- 
emitting neutron stars and currently a few tens of these ob- 
jects are known ( see th e RRATalog]). In the discovery paper, 
iMcLaughlin et alj ( I2OO6I) attempt to simulate the size of the pop- 
ulation based on the 1 1 sources they report. Nevertheless, no syn- 
thesis models have been made to describe in detail their individual 
characteristics since our current understanding of the population as 
a whole is uncertain. To build our model, we use the measured pe- 
riods and other parameters of more than 50 known RRATs reported 



^ http://www.atnf.csiro.au/research/pulsaiypsrcat/ 
* http://www.as.wvu.edu/~pulsar/rratalog/ 



by 'McLaughlin et al.' ('2006'), 'Hessels et al.' (2008'), 'Oeneva et alj 
(2009), Keane et al. (2010) and Burke-Spolaor & Bailes (201(F 
With this empirical model, we aim to describe the general prop- 
erties of a galactic population of about 40,000 RRATs, in order 
to estimate the distribution of periods N{P), luminosities L(P) and 
pulse rate N. For the period distribution we use the following ana- 
lytic function: 



N{P) = A \og{P) + B, 



(6) 



while for the luminosity we used the following relation: 

logL = Dlog(P) + £ + r. (7) 

with r being a normal distribution with zero mean and standard 
deviation of 0.4 (obtained from our best fits). The sporadic nature 
of RRATs can be quantified by the pulse rate N of the individual 
objects. This parameter was modeled according to: 



log/V = Flog(/') + G + Y 



(8) 



In the last equation Y is a normal distribution with zero mean and 
standard deviation 0.7 (as above, obtained from best fits). The fitted 
parameters are tabulated in the fifth column of Table |4] For each 
RRAT in our model we randomly assigned a distribution for the 
intensity of its individual pulses with the same three types of pulse- 
energy distributions used for pulsars. We have scaled the luminosi- 
ties in the same way as described in §5.1.1. The spectral index of 
RRATs is still poorly known, therefore we have simulated popula- 
tions assuming an average spectral index like pulsars (a = -1.82). 
We note that all these parameters are uncertain, given the small 
number of known sources. 



5.7..? Monte-Carlo Simulations 

Any pulsar located in M31 must be very luminous in order to be 
detected. Thus, a search for pulsars in a whole galaxy such as M31 
is an opportunity to set upper limits on the maximum luminosity of 
radio pulsars in general. To quantify this luminosity, we performed 
different simulations of both pulsars and RRATs. Each simulation 
had different luminosity distributions, covering a range from that 
observed in our Galaxy to a population 10 times brighter, in order to 
establish which luminosity distributions would ensure a detection 
with the WSRT From now on, a model realization is defined as a 
single population of pulsars. For each simulation, we produced 100 
pulsar/RRAT populations. For each of these models, and based on 
the luminosity distributions, we calculated the mean flux of each 
pulsar given the distance to M31. We then compared these mean 
fluxes with the sensitivity of our periodicity searches. Similarly we 
have simulated distributions of individual pulses for each pulsar 
and RRAT. The brightest pulses of each distribution were compared 
with our single pulse search sensitivity. 

A detection in our models is any pulsar or RRAT that has a 
flux above the minimum thresholds given in §3.6. Those pulsars 
or RRATs that were above the minimum threshold were recorded 
along with their fluxes S and pulse-energy distributions. The re- 
sults are summarized in Tables [5] and [6] respectively. They show in 
each column the average number of simulated pulsars detected for 
each assumed distribution. In the case of single-pulse detections, 
we also show the details of the power law, normal, or log-normal 
distributions for those cases where we had a detection. A general 
result of our simulations is that normal radio pulsars can be more 
easily detected than RRAT-like objects. 

Table [5] shows that we should be able to detect at least a few 
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Figure 10. One example model realization for radio pulsars. The left panel shows a histogram of the number of pulsars versus peiiod. Our models cover periods 
ranging from a few tens of milliseconds to a few seconds. The right panel shows the distribution of luminosities versus period. The dotted line represents the 
most luminous known pulsar in our Galaxy, PSR B 1 302-64, which has a L328 ~ 37 Jy kpc- . 
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Figure 11. One example model realization for RRATs. As in Fig. 1101 the left panel shows the distribution of the number of RRATs versus period, based on the 
known periods of RRATs. The right panel shows the distribution of luminosities versus period. The dotted line shows the most luminous RRAT known in our 
Galaxy, J1819-1458, with peak luminosity L328 ~ 5170 mjy kpc". The luminosities shown here were scaled with frequency as explained in the text. 



Table 4. Summary of the fitted parameters from our pulsar and RRAT models. The various parameters used are described in Eqs. 4-8. The columns marked 
with d^^ are the model fits using the standard inverse square law for flux as a function of distance. The columns marked with d^^ and d^'-^ show the fits for 
other assumed flux scaling laws, as explained in §5.2. The spaces marked with "-" indicate that the models did not require those parameters. 



Parameter 




Pulsar 






RRATs 




A 


0.510 


0.510 


0.510 


-0.467141 


-0.467141 


-0.467141 


B 


2.710 


2.710 


2.710 


679.3115 


679.3115 


679.3115 


C 


0.340 


0.340 


0.340 








D 


-0.230 


-0.211798 


-0.1931 


-0.045 


-0.0742648 


-0.0645126 


E 


3.212 


2.36363 


-0.0645 


2.590 


3.47843 


3.18352 


F 








-0.2633 


-0.2633 


-0.2633 


G 








1.91755 


1.91755 


1.91755 
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Table 5. Monte-Carlo simulations for different luminosity distribution scalings (shown in the first column) for pulsars. In the other columns we show the 
number of pulsars detected using the different methods, namely periodicity and single pulse searches. We did simulations using shallow (o- ~ 1) and steep 
{a > 2) power law indices. For each case we show the number of detectable pulsars using a flux law that scales as d"^, rf"'-^ and t/"' (see text in §5.2). The 
spaces marked with "-" indicate that we did not model those conditions since they were intermediate between the lowest and highest luminosity models. 
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lOxL 5 375 14565 222646050 



Table 6. Monte-Carlo simulations made for RRATs using different luminosity laws (first column). L is the luminosity as seen for the population of RRATs as 
observed in our Galaxy (L = Sd^). In the second and third columns we show the number of RRATs detected via their single pulses. We also show the number 
of detections for each of the flux distribution models. In the power law models, we did simulations using shallow (a ~ 1) and steep (a > 2) power law indices. 
As in Table [s] the numbers in each column below the distributions show the expected number of detectable RRATs using a flux law that scales as d^^, d^^ ^ 
and d"' (see text in §.5.2). 
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single pulses from normal radio pulsars located at the distance of 
M31 only if the intensity distribution of their individual pulses 
shows either a shallow power law (1< o- ^2) or a log-normal dis- 
tribution. This is the case for all known pulsars, which have fluxes 
that scale proportional to cT^. As can also be seen in Table [5] if 
the population of pulsars in M31 would be at least an average of 
5 times more luminous than the one in our Galaxy, our periodic- 
ity searches might detect a similar number of pulsars as our sin- 
gle pulse searches. Thus PSR B 1302-64 with L j,2g = 37 Jy kpc^, 
the m ost luminous known pulsar in our Galaxy ^Manchester et al.l 
Il978l) , could be detected in M31 if it were 5 times brighter. Based 
on our model of the luminosity distribution, RRATs located in M3 1 
are less likely to be found than pulsars, and only those populations 
exhibiting more than 5 times the observed luminosity of RRATs 
in our Galaxy, and having individual pulses with fluxes distributed 
with either a log-normal or a shallow power law could be detected. 
We note that this result was obtained by assuming that the aver- 
age spectral index of RRATs is the same as that for normal radio 
pulsars. If the average index is a factor of 2 less steep we might 
detect at least 1 of these objects in M31. We note that a source like 
J1819-1458, the m ost luminous R RAT known with L328 = 5170 
mjy kpc^ ( McLaughl in et al .I2OO6I) would have to emit pulses about 
25 times brighter to be detected in M3 1 . It is important to mention 
here that the dwell time of the observations plays a key role in the 
detection of a single pulse; the longer we observe the higher the 
chances of observing a bright pulse. 



The fact that we have not detected any pulsars in our periodic- 
ity searches can be attributed to a few reasons. The first could sim- 
ply be that pulsars are not luminous enough to be detected at that 
distance with the WSRT. If that is the case, these observations im- 
ply that there are no pulsars with luminosities L328 > 185 Jy kpc^. 
Secondly, the population of pulsars in M3 1 could be smaller than 
that of our Galaxy. This argument is supported by observations 
made at other wavelengths which show that M3 1 has a star forma- 
tion rate of 0.35 M© yr"' ( Wi!li_ams 2003 and lBarmbv et al.ll2006l) 
about a factor of 3 lower than in our Galaxy. Even though the num- 
bers quoted here are uncertain, this suggests that the number of 
neutron stars recently produced in M31 is less than in our Galaxy. 
Further evidence in favor of a smaller population of radio pulsars 
comes from the low me tallicity observed in M3 1 , a factor of 2 lowe r 
than that in our Galaxy iSmartt et al.l200ll and lTrundle et alj2002h . 
Photometric surveys made in the disk of M3 1 have revealed a large 
population of old metal-poor stars (Renda et al. 2005) which may 
mainly form bla ck holes when they collapse as supernova (see 
iHeger et al. I [2003) . Therefore the population of pulsars in M31 
might also be older than in our Galaxy, with most of its members 
being too faint to be detected. Beaming could also cause us to miss 
the few pulsar candidates found bv iGelfand et al.i t2005i) . 
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5.2 Constraining the d ^ flux law for radio pulsars 

The "pseudo-luminosity" of a pulsar at a given observing frequency 
(Ly) and the distance to the pulsar (d in kpc) relate as: 5y oc Ly/d-. 
The mean flux density 5y is based on the integrated intensity un- 
der the pulse, averaged over the full period. The most widely used 
method for determi ning distance is the "D M-distance", based on 
the NE2001 model dCordes & Lazidl2002 ) for the distribution of 
free electrons in our Galaxy (ne). This model allows distances to be 
estimated to an accuracy of about 20% on average for large sam- 
ples but can give an error of 50% or more for individual objects. 
Since the model describes the entire distribution of free electrons 
in our Galaxy, it is used to estimate distances to pulsars from their 

mea sured DM. 

ISingleton et alj j2009l) made a maximum likelihood analysis 
using the observed fluxes and DM-distances of pulsars. The con- 
clusions of their analysis are that either i) the distances are radically 
and consistently incorrect by factors of 10, or ii) pulsars have a flux 
that falls off more slowly with distance than expected, following 
a oc d^^ or oc t/"'-^ power law. This discrepancy in fluxes is ex- 
plained by the presence of a component in the flux du e to the emis- 
sion of superluminal polarization currents ( lArdavan~e t al. 2008). In 
order to test this claim, we performed simulations like those de- 
scribed in §5.1, but assuming that the flux scales as Sy oc Ly/d^ 
with e = 1, e = 1.5 and e = 2, in order to predict the number 
of detectable pulsars at the distance of M31. Our detailed simula- 
tions paid particular attention to the luminosity of the sources and 
a range of different models for the single pulse intensity distribu- 
tions. The results are summarized in the Tables[5]and|6]for pulsars 
and RRATs respectively. Our results clearly show that if the flux 
density obeys a d"' law, then we should have detected many hun- 
dreds of pulsars. The fact that our observations and those made of 
M33 jMcLaughlin & Cordesir2003h revealed only a handful of sin- 
gl e bright pulses, and th at our simulations using the fluxes scaled 
as ISingleton et all l l2009l) show a large number of detections, lead 
us to conclude that the flux of pulsars does not scale as d^' or d^' ^ 
but simply as d^^ as always assumed (unless the pulsars in M31 
and M33 are systematically far weaker than in our Galaxy). This 
argues that some other effect than distance errors o r a modified lu- 
minos ity law is responsible for the effect seen bv ISingleton et al.l 
(12009). We suggest that the observed deviation may be caused by 
an unaccounted selection eff^ect in the flux measurements of radio 
pulsar surveys. 

5.3 What is the detection at DM=54.7 pc cm^? 

The dispersed nature of the pulses detected at DM=54.7 pc cm"' 
and their similar properties (intrinsic width and peak flux) is con- 
sistent with individual bright pulses observed in known pulsars. 
The individual pulses show peak flux densities between 2.8-3.2 
Jy and a luminosity > 30 Jy kpc^ assuming they are located in 
our Galaxy at the distance predicted by the NE2001 of d = 2.9^2 1 
kpc. However given that the DM contribution of our Galaxy may 
be as low as 45 pc cm"' in this direction, this makes it possible that 
these pulses originate in M31, in which case these pulses would be 
~ 70000X more luminous. Neither the arrival-time distribution of 
individual pulses shown by the candidate at DM=54.7 pc cm"' or 
the intrinsic width of the pulses is consist ent with that observed 
in the Crab 's giant pulses, as reporte d in |m cLaughlin & Corded 
( l2003h and iKaruppusamv et al] ( l2010t) . If these pulses would be 
giant pulses from a distant Crab-like pulsar, then we would have 
detected a similar or higher number of pulses per hour (N> 10 h"' 



for the Crab) in many, if not all, of our observations. Assuming the 
signals are indeed astronomical in nature, it is thus more likely that 
the source is showing normal pulsar or RRAT-like pulses in a burst, 
as has been seen in other, known sources. 

5.4 Detectability of cosmological millisecond radio bursts 

Though we did not obviously detect any such events, the long 
dwell times of our observations provided the potential to also de- 
tect extragalactic millisecond radio bursts occurring at cosmologi- 
cal di stances and in the general direction of M31 (cf. lLorimer et al.l 
|2007|) . Here we estimate the number of such bursts that could be 
observable in our data. Lorimer et al. (2007) estimated the rate of 
such events by using the properties of the Parkes Multibeam Sur- 
vey over an area of 5 square degrees (1/8250 of the entire sky) at 
any given time over a 20-day period. Assuming the bursts to be 
distributed isotropically over the sky, they infer a nominal rate of 
8250/20 ~ 225 similar events per day. Based on those estimates, 
the expected number of bright single bursts per day per square de- 
gree, «B, is ~ 0.01 bursts day"' sqr deg"'. 

We have observed for a total of 80 h (3.33 days) with a FoV 
9 = 2.57° across. The area of the sky covered by the primary beam 
is thus Q = k[l - cos(0/2)] = 5.15 sqr. deg., with k = 2- 3283 • n. 
Therefore, we estimate the total expected number of pulses to be: 

0.01 

= X 3.33 day x 5.15 sqr. deg. = 0.17. (9) 

day • sqr. deg. 

This implies that we would likely need to observe longer to de- 
tect such an extragalactic millisecond radio burst, assuming that 
a cosmological radio burst can still be detected at 328 MHz af- 
ter crossing near M31. The bright single pulses listed in Table 
All are arguably not like the millisecond radio burst detected by 
Lorimer et al. |(|2003), jecause i) they are not nearly as bright as the 
one reported there, ii) we do not have precise localization of these 
sources and hence they may well originate in M31, and iii) they 
are likely not at cosmological distances due to the DMs observed, 
which can be easily explained by the galactic medium of M3 1 . 



CONCLUSIONS 

We have presented the most sensitive search for pulsars and radio 
transients ever made in the direction of M3 1 . We present our con- 
clusions below. 

(i) We have detected a handful of single bright pulses with a DM 
that locates them beyond our Galaxy and plausibly within M31. 
Moreover we detected a few repeating bright pulses that consis- 
tently show the same pulse width and DM. Some of them have a 
DM consistent with that of a source at the distance of M3 1 while 
others could be located within our own Galaxy. 

(ii) We performed simulations using a simple synthesis code 
to estimate the properties of a population of pulsars and RRATs 
in M31. From our simulations, we conclude that the only means 
to detect pulsars in M31 with the sensitivity of the WSRT in the 
8gr8 mode is by detecting bright pulses described by a shallow 
power law or a log-normal distribution and assuming that the pop- 
ulation of pulsars in M3 1 has similar characteristics to that of our 
Galaxy. The detections reported here might be associated with pul- 
sars or RRATs that show these characteristics. Normal pulsars are 
more easily detected than RRATs because the individual pulses of 
RRATs are generally no brighter than normal pulsar pulses. 
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(iii) As a consequence of observing with the WSRT in the 8gr8 
mode, we were not able to locate sporadic bursts to high accuracy 
within the primary beam. Further observations with shorter dwell 
times made with the GBT did not confirm these detections at the 
expected locations. 

(iv) The fact that we do not detect any pulsars in our periodic- 
ity searches can be explained by the absence of extremely bright 
pulsars (L328 > 185 Jy kpc^) and/or a smaller population of pulsars 
and RRATs in M3 1 compared with our Galaxy. 

(v) Our observations indicate that t he flux of pulsars doe s not 
scale as rf"' or J"'^ (as suggested by ISingleton et al]|2009l) . but 
simply as d"^. Nevertheless there is an inconsistency that pos- 
sibly can be attributed to selection effects in the pulsar surveys 
which show a bias towards the detection of pulses with low fluxes 
(Malmquist bias) or perhaps a combination of this effect with the 
distance of pulsars. 

In the future, WSRT, in combination with the APERTIF fo- 
cal plane array, will allow for similar surveys to the one presented 
here, but with a larger instantaneous FoV and a higher observing 
frequency (1.4 GHz, which is likely better for detecting highly 
dispersed pulses). Future searc hes of M31 with LOFAR are also 
promising jStappers et al.l20lll) . and the SKA is likely to make the 
first detections of pulsars in the nearby galaxies visible from the 
Southern Hemisphere. 
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Table Al. Single bright pulses observed in the direction of M3 1 with S/N> 7.0. The pulses shown in this table were only seen once 
at this spatial location and DM. We indicate here the number of sub-beams in which each pulse was detected, the peak DM of each 
detection, the distance (using the NE2001 model), S/N, approximate pulse width (W), the arrival times relative to the start of the 
observation and the start time of each observation (topocentric). The last column shows the name and the number of the pointing in 
which the detections occurred (for example PNTOl 3 is the third pointing out of four pointings made with the name PNTOl etc). The 
detections in boldface are shown in Figs. lA2l and lAll and are representative examples of pulses with DM > 100 pc cm"'. 
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Figure Al. Examples of the events detected as single bright pulses, which are summarized in Table IaT] The panels show from top to bottom single bright 
pulses detected at DM = 104 pc cm"^ (t ~ 3436 s) DM =264 pc cm"' (t ~ 2506 s) and DM = 267 pc cm"' (t ~ 2541 s). In the left panels we show, for clarity, 
the peak DM of each detection with horizontal dashed lines. The high DM of these detections suggest that they are located outside our Galaxy. The right panels 
show the S/N vs. DM plots. The x axes of the right panels go only to t = 3500 s. 
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Figure A2. Pulse profile of a single pulse detected with a DM= 139.41 pc cm"^. The upper panel shows the time series at the DM of the detection. The lower 
panels show how the individual pulses decrease in S/N at adjacent DM trials both above and below DM=139.4 pc cm"^ (thick line). The high dispersion of 
this pulse locates it at a distance > 48 kpc, therefore the pulse plausibly originates in M31. The details of this detection are shown in Table lAll 

Table Bl. Bright pulses detected more than once at the same DM during 80 h of observations of M31. We show the same parameters 
as in Table lAll including the anival times of the detections (topocentric). We stress that these candidates are potentially detected from 
the same sky direction. 
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Figure Bl. Pulse profile of two pulses detected at DM=72.4 pc cm"' at dilferent times, as shown in Table Bl. For the two pulses, the upper panels show 
the time series at the DM of the detections. The lower panels show how the individual pulses decrease in S/N at adjacent DM trials both above and below 
DM=72.4 pc cm~' (thick line). The high DMs of these pulses locates them at a distance > 48 kpc; therefore these pulses plausibly originate in M31. The 
details of this detection are shown in boldface in Table lBll 
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